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An electrical contact resistance (ECR) theory is presented for conductive and rough (fractal) 
surfaces separated by a thin insulating film, which is treated as an energy barrier that impedes 
current flow due to the electric-tunnel effect. The analysis yields insight into the effects of film 
properties, current flow, surface topography, mechanical properties, and contact load on the ECR. It 
is shown that the variation of ECR with the contact load is less pronounced than that observed in the 
absence of an insulating layer due to the intrinsic voltage dependence of the tunnel resistance and 
the associated voltage compensation mechanism. The effect of nonohmic behavior on the 
relationships of the ECR with the contact load and the real contact area is discussed and results are 
compared with approximate analytical relationships developed herein. The relationship between the 
real contact area and the ECR depends on the current intensity and film properties and is 
independent of the surface topography and mechanical properties. Approaches for determining the 
surface roughness, mechanical properties, insulating film properties, and real contact area from ECR 
measurements are interpreted in light of the theory developed. © 2004 American Institute of 


Physics. [DOI: 10.1063/1.1629392] 


I. INTRODUCTION 


Thin insulating films formed at the contact interfaces of 
components that consist of conductive materials play an im- 
portant role in the performance of various electromechanical 
devices,!~? such as circuit breakers, connectors, relays, and 
switches. Significant advances in the development of reliable 
microrelays* and microswitches> that can be actuated elec- 
trostatically have been recorded recently. Although the con- 
tacting surfaces of these microdevices consist of structural 
polysilicon layers metallized with gold, there is a growing 
concern about the presence of contaminants and insulating 
films at contact interfaces that may greatly increase the elec- 
trical contact resistance (ECR). For example, the formation 
of an insulating oxide layer on palladium, which is often 
used in small-size relays, has been reported to be the main 
cause of electrical contact deterioration.° Thus, there is a 
great demand for new materials that can increase the device 
lifetime, while maintaining a low and stable contact resis- 
tance for a long period.’ 

The high tunnel resistance due to an insulating thin film 
may degrade the performance of components produced with 
various bonding techniques, such as conductive adhesive 
joints,” flip-chip bonding and ball-grid array sockets,!° and 
resistance spot welding.'! Moreover, the existence of native 
insulating oxides prevents the use of certain materials in 
joining applications. ! Despite the fact that in most applica- 
tions low ECR is desirable for better performance,'? there are 
also situations where high ECR enhances the device 
lifetime.'* The strong effect of an insulating film on ECR is 
beneficial in many applications, such as hydrogen 
detection," determination of the properties of surface 
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oxides, !® and detection of film formation at sliding 
interfaces." Therefore, it is important to understand the role 
of a thin insulating film in ECR measurements. 

One of the early ECR models that takes into account the 
presence of a thin film adsorbed at the surfaces is attributed 
to Holm, ! who observed that the ECR is affected by both 
the constriction resistance and the film resistance. A general- 
ized ECR theory!” that uses fractal geometry for a descrip- 
tion of the surface topography,” elastic-plastic deformation 
of the contacting asperities,”! and size-dependent electrical 
constriction resistance of microcontacts”” comprising the real 
contact area was developed for homogeneous conductive 
surfaces. In general, theoretical predictions for homogeneous 
conductive surfaces underestimate the ECR measured in the 
presence of thin insulating films at the surfaces.” In a pre- 
vious ECR study,” a surface thin film was modeled as an 
ideal insulating material. However, when two surfaces are 
separated by a sufficiently thin insulating film, it is possible 
for current to flow due to the tunnel effect.” Another ECR 
study that accounts for the presence of interfacial insulating 
thin layers”° was based on the general theory of the electric- 
tunnel effect between similar electrodes separated by a thin 
insulating film.” An equivalent electric circuit for the con- 
tact interface was proposed”° and the current intensity and 
voltage drop across the interface were found to be in good 
agreement with experimental results. However, although the 
equivalent electrical circuit proposed in Ref. 26 characterizes 
the global behavior at the contact interface, it does not ac- 
count for the physical behavior of the microcontacts. 

Another shortcoming of a previous ECR model” that 
was based on detailed analysis at the asperity microcontact 
level is the use of a statistical approach to represent the sur- 
face topography.”° This approach is limited by the depen- 
dence of statistical roughness parameters on the sampling 


© 2004 American Institute of Physics 


OL-2r:vL EZOZ Jequisceq €; 


J. Appl. Phys., Vol. 95, No. 2, 15 January 2004 


length and the resolution of the measuring instrument and, 
therefore, cannot yield unbiased information about the sur- 
face topography, which plays an important role in ECR mea- 
surements. This deficiency can be overcome by using a frac- 
tal topography description that is based on scale-invariant 
parameters.” Fractal geometry provides an effective means 
by which to model engineering surfaces that exhibit random, 
multiscale topographies due to its intrinsic advantages of 
scale invariance and self-affinity.*” A fractal contact analysis 
that accounts for the elastic, elastic-plastic, and fully plastic 
deformation of the asperity microcontacts at rough surfaces 
that exhibit fractal behavior was derived in a previous 
study.7! 

Although earlier ECR studies have yielded valuable in- 
sight into the origins of ECR, a general theory for conduc- 
tive, elastic-plastic, and rough surfaces separated by a thin 
insulating film has not been developed yet. Thus, the main 
objective of this study was to produce a general ECR theory 
for conductive surfaces separated by a thin insulating film. 
The present ECR theory is based on fractal geometry to de- 
scribe the surface topography,” elastic, elastic-plastic, and 
fully plastic deformation of asperity microcontacts,”! and 
quantum mechanics considerations for the electric-tunnel ef- 
fect through the thin insulating film.” It is demonstrated that 
the constriction resistance is negligible with respect to the 
tunnel resistance. This means that the ECR behavior can be 
profoundly changed in the presence of an insulating thin film 
at the contact interface. Consequently, there is a marked dif- 
ference with the theory derived for conductive rough 
surfaces, |? where the constriction resistance is the only 
mechanism that impairs current flow. 


Il. THEORETICAL TREATMENT 
A. Description of the surface topography 


Fractal geometry”? provides a means by which to over- 
come limitations associated with the resolution of the instru- 
ment used to measure surface roughness parameters. Two 
important parameters in a description of fractal topography 
are the fractal roughness G and the fractal dimension D 
(2<D<3). The fractal roughness G is a height scaling pa- 
rameter that is independent of the frequency. A rougher sur- 
face is characterized by higher G values. The fractal dimen- 
sion D determines the contribution of high- and low- 
frequency components in the surface profile. Hence, high 
values of D indicate that high-frequency components are 
more dominant in the surface profile than low-frequency 
components. For fixed fractal roughness G, higher values of 
the fractal dimension D yield smoother topographies; how- 
ever, the amplitude ratio of high-to-low frequencies of the 
surface increases with D. Details about characterization of 
the surface topography using fractal geometry can be found 
elsewhere? -2129-31 and, therefore, are not repeated here for 
brevity. 


B. Mechanical contact model 


Contact of two rough surfaces [Fig. 1(a)] is equivalent to 
contact of a smooth (flat) half-space with reduced elastic 
modulus E=[(1— v?)/E1+(1-—1Ż)/E,] !, where v, and 
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(b) 
| 
(c) 


FIG. 1. (a) Rough surfaces in normal contact separated by an insulating thin 
film, (b) equivalent contact model, and (c) electrical analog of contact resis- 
tance. 


v, and E; and E, are the Poisson ratios and elastic moduli of 
the two surfaces, respectively, and a rigid rough surface with 
topography equivalent to those of the two rough surfaces* 
[Fig. 1(b)]. For fractal surfaces, this implies that the power 
spectrum of the equivalent rough surface is equal to the sum 
of the power spectra of the two surfaces.*! The equivalent 
thickness ¢ of an insulating film covering the surfaces is 
t=t,+t2, where tf, and t, denote the thickness of the film on 
each surface. For very thin films (i.e., t<50 Aj," in which 
the electric-tunnel effect occurs, the contact mechanics of 
layered media are dominated by the mechanical properties of 
the substrate.”! Surface contact produces numerous circular 
asperity microcontacts [Fig. 1(b)], assumed to be sufficiently 
separated from each other in order for asperity interactions to 
be secondary. This assumption is reasonable for lightly 
loaded contacts where the real contact area A is a small frac- 
tion of the apparent contact area A, .°8 Based on these as- 
sumptions and on the calculated mean contact pressure and 
contact area of the asperity microcontacts, the total contact 
load P and real contact area can be obtained using an inte- 
gration procedure that accounts for microcontact contribu- 
tions to the previous contact parameters.7°7! 

For the case of fully plastic deformation of all the mi- 
crocontacts, an explicit relationship between the real contact 
area and the contact load can be obtained because the mean 
contact pressure is equal to the hardness of the softer surface 
H, assumed equal to cY, where Y is the yield strength and c 
is typically equal to 2.8.3 Hence, the following relationship 
that is independent of surface topography is obtained:*4 


A*=(E/H)P*, (1) 
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where the dimensionless real contact area and contact load 
are defined as A*=A/A, and P*=P/(A,E), respectively. 

For elastic or fully plastic microcontact deformation, the 
total contact load and real contact area, determined analyti- 
cally by integrating the contribution of the elastic and fully 
plastic microcontacts, are given by” 


.\ (3—D)2 (3-D)/2 
= (75 a aç” 


I% I% 
ay ay, 


(11-—2D)/2 _ 
2 | D-1 Jos y)?GrO-2) 


3 q(4-D2\5—2D 
r% \ (5-2D)/2 
xaea] e, O 
a,” 
and 
qt D-1 
~ 2(3-D) 
rx \ (3—D)/2 rx \ (3-D)/2 
xai*| | =e -2| 25 +1|. 0) 
ar, t$ Ik 7 
L ay, 


The dimensionless fractal roughness G* is defined as G* 
=G/A va , yis a scaling parameter used in fractal description 
of the surface topography (typically, y=1.5),”° and a% and a} 
are the smallest and largest truncated microcontact areas, 
respectively”! The critical truncated microcontact area a A 
separating the elastic from the fully plastic deformation re- 
gime is given by 


1(D—2) 


> 


2 
a= 29 2D) (P Db 1 GP (5 Iny 


(4) 


where b=[77(0.454+0.41,)/2]’, in which v; is the Pois- 
son ratio of the softer material. Microcontacts with truncated 
contact areas a’>a/, and a’ <a’, are in the elastic and fully 
plastic deformation regimes, respectively.”” The above mi- 
crocontact areas are normalized with respect to A, and the 
dimensionless parameters are denoted by an asterisk. 
Relationships for the mean contact pressure and contact 
area have been derived for a single asperity in the elastic, 
elastic-plastic, and fully plastic deformation regimes.”! How- 
ever, because the constitutive relationships for elastic-plastic 
deformation are more complex than those of the simpler 
elastic and fully plastic deformation,” a numerical integra- 
tion scheme was developed in order to calculate the total 
contact load and real contact area.” Consequently, the total 
contact load and real contact area can be expressed as 


N(ag) 
P=) _ F,, (5) 
i=l 
and 
N(ag) 
A= 3 ai, (6) 


i=1 
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where N(aç) is the number of truncated asperities that pos- 
sess areas greater than the smallest truncated contact area 
ay. The contact load F; and contact area a; of the ith micro- 
contact are calculated using the appropriate constitutive 
relationships.”! 


C. Electrical contact resistance at a single 
microcontact 


The ECR at each microcontact R; consists of the con- 
striction resistance R.,; due to the convergence and diver- 
gence of current flow'® and the tunnel resistance R si due to 
the presence of an insulating film that introduces a potential 
barrier that impedes the flow of electrons.” An electrical 
analog of the total contact resistance due to the constriction 
resistance produced by the microcontacts comprising the real 
contact area and the insulating film between the contacting 
surfaces is presented in Fig. l(c). It has been shown! that, 
for lightly loaded surfaces, where the radius of each micro- 
contact r; is smaller than the average electron mean-free path 
of the contacting surfaces, \=(A,;+),)/2, the constriction 
resistance is dominated by the Sharvin mechanism. Accord- 
ing to this mechanism, the electrons travel through the mi- 
crocontacts without undergoing any scattering, and the con- 
striction resistance is given by” 


Ram- (7) 


where p is the average specific resistivity of the contacting 
surfaces, p= (pı + p2)/2. 

For electrons to pass from one surface to the other in the 
presence of an insulating film they must have sufficient en- 
ergy to surmount the barrier produced by the insulating film, 
i.e., they must enter the conduction band of the insulator, a 
process known as thermionic emission. According to classi- 
cal physics, the electrons cannot penetrate the barrier if the 
electron energy is less than the height of the interfacial bar- 
rier. However, according to quantum theory, a finite probabil- 
ity exists for the electrons to tunnel through the barrier, de- 
pending on the size and shape of the barrier encountered by 
the electrons. 

Simmons?” developed a formula for electric tunneling 
through a potential barrier of arbitrary shape in a thin insu- 
lating film that separates two similar electrodes. The formula 
was derived for a rectangular barrier and included the image 
force. The effect of the image force is to reduce the area of 
the potential barrier by rounding off the corners of its distri- 
bution, thereby reducing the thickness of the barrier and, 
consequently, increasing the current flow between the elec- 
trodes. The analysis presented in Ref. 27 is for low tempera- 
tures, where the thermal current can be neglected, thus re- 
stricting electron transport between electrodes by the tunnel 
phenomenon. However, it has been shown! that the equa- 
tions derived in Ref. 27 are also applicable at higher tem- 
peratures. Consider a single microcontact of area a; [Fig. 
1(b)] covered by a thin insulating film of thickness t, dielec- 
tric constant K, and energy height above the Fermi level of 
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the conductive surfaces œo. The current intensity 7; through 
the microcontact and associated voltage drop V; in three dif- 
ferent voltage regimes are given by?” 


1,=(3.16X 10!/AS) o! 


<exp(—1.025AS@;")V,a;, for V;=0, (8) 
where 
Sa(t—S1) 

L= go- [5.75/(KAS) In h, (9) 
and 

AS=S,-S], (10a) 
where Sı and S, are given by 

S,=6/(K go), (10b) 
and 

S =t—6/(K@o). (10c) 


Tunnel resistance R,; for V;=0 is derived from Eq. (8) to be 


V; AS exp(1.025AS@;”) 1 
OL, 3.16X10% pl? a; 


(11) 


Equation (11) shows that R,; is independent of the voltage or 
the current applied to the microcontact to measure the tunnel 
resistance and, therefore, is ohmic resistance. As can be seen 
from Eqs. (9)—(11), R,; is inversely proportional to a; and 
increases with t, K, and Pp." 

For a drop in voltage across a microcontact V; that is not 
very small,” 


1;=[6.2X 10!°/(AS)?]{ @, exp(— 1.025AS 97”) 
—(p_tV;)exp[— 1.025AS(9,+ V) ]}a;, (12) 


where 
PL= Go-LVi/(2t) (Si +S2) 
—[5.75( KAS) ]In ae, (13) 
Si (t—S>) 
and 
S,=6(K go), (14) 
S2=t[1 —46/(3 poKt+20—2V;Kt)]+6/(K po), 
for V;<@o, (15a) 
or 
S,=(@oKt—28)/KV;, for V;>@o. (15b) 


As can be seen from Eq. (12) the relationship between cur- 
rent and voltage is nonlinear. Therefore, it is not possible to 
derive an explicit relationship for the tunnel resistance, such 
as that given by Eq. (11). Moreover, because z; and AS are 
functions of the voltage, the tunnel resistance is nonohmic 
and decreases with an increase in voltage.” 

In Eqs. (8)—(15), and throughout this study, J; is ex- 
pressed in A, V; and @p in V, t, S4, and S, in A, and a; (as 
well as all other areas) in cm’, 
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Ill. APROXIMATE ANALYSIS 


To obtain insight into the contact electromechanics be- 
havior, it is instructive to first consider the following analysis 
derived for the low-voltage regime (V;=0) and either elastic 
and fully plastic microcontacts or solely fully plastic micro- 
contacts. Figure l(c) shows that the constriction resistance 
and the tunnel resistance of a microcontact are in series. The 
ratio between the restriction resistance and the tunnel resis- 
tance in the low-voltage regime (V;=0) is obtained by di- 
viding Eq. (7) by Eq. (11), 


Ra 3-16 10° @7"Ap 
Ri AS exp(1.025AS 9”) , 


(16) 


When R.;/R,;<1, the constriction resistance can be ne- 
glected in favor of the tunnel resistance. The ratio R,;/R;; is 
examined for a conductive surface that possesses a fairly 
high resistivity of p= 1074 Q cm and \=3 X10 °cm, and a 
relatively low energy barrier with t= 10 A, o=1 V, and 
K=4. For this case, Eqs. (9), (10), and (16) yield R,,/R,; 
=2.7X1073, indicating that the constriction resistance can 
indeed be neglected as secondary; hence, the contact resis- 
tance at a microcontact is R;=R,;. This is a reasonable ap- 
proximation for conductive surfaces separated by a thin in- 
sulating film, especially if the surfaces are very conductive 
and the energy barrier imposed by the insulating film is 
higher and/or wider. Equation (16) is valid when V;=0, i.e., 
when the tunnel resistance is ohmic. For higher voltages, the 
tunnel resistance [Eqs. (12)—(15)] decreases with an increas- 
ing drop in voltage? and, hence, the ratio R,.;/R,; increases. 
As shown in Ref. 27, the tunnel resistance decreases by ap- 
proximately an order of magnitude when the voltage in- 
creases from 0 to 1 V. However, despite this effect, the ratio 
R,,/R,; assumes values significantly less than 1. 

As mentioned in Sec. IIB, the equivalent contact model 
of two rough surfaces consists of a smooth elastic-plastic 
medium in contact with a rigid rough surface. The corre- 
sponding real contact area is the sum of the discrete asperity 
microcontacts. The total ECR is assumed to be the sum of 
individual parallel resistances corresponding to the restric- 
tion resistances of individual microcontacts. 

The electrical conductivity C; of a single microcontact in 
the low-voltage regime (V;=0) is obtained from Eq. (11) to 
be 


a 3.16 100g]? exp(— 1.025AS¢7”)a; 
i i AS i 


Using an integration procedure for elastic and fully plastic 
19 


(17) 


microcontacts, ~ the total electrical conductivity C is ob- 
tained as 
a, ar 
g=] Cipla’ )n(a' aa’ + | C;-(a')n(a')da', (18) 
ag al 


where n(a’) is the size distribution function of truncated 
microcontacts,”” and subscripts p and e denote plastic and 
elastic, respectively. The electrical conductivity of a single 
microcontact in the elastic and fully plastic deformation re- 
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gimes, C;, and C;,, respectively, can be obtained by substi- 
tuting the appropriate contact area relationships”? into Eq. 
(17). After integration, Eq. (18) yields 


ee 3.16X 10!" exp(—1.025AS@;") D-1 
AS 2(3-D) 


1\ (3-D)2 ,\ (3-D)/2 
a as 
— —2|— +1]. 
ay ay 


The dimensionless electrical conductivity C* is defined as 


Xar (19) 


s Ct 
3.16X 100A, 
1/2 
_|ĦL t 1/2 _ 
(2 PE 1.025AS 9; )53=D) 
(3—D)/2 (3-D)/2 
al* al* 
x<al* | =) =, a +1), (20) 
ar ar 
and, thus, the dimensionless ECR can be expressed as 
3.16X 10V oA, 
R*=(C*)"'= R. (21) 


t 


From Eqs. (4), (9), (10), (20), and (21) it can be seen that 
the dimensionless ECR depends on fractal parameters D and 
G*, mechanical properties E/Y and vı, insulating film 
thickness t and properties K and go, and dimensionless 
smallest and largest truncated microcontact areas aç* and 
a; , respectively. For a continuum description, ag must be 
greater than the atomic dimensions, e.g., about six times the 
lattice dimension.”! For given surface roughness, mechanical 
properties, and insulating film properties, the only unknown 
parameter in Eq. (20) is a;*, which can be found implicitly 
from Eq. (2) as a function of contact load P*. Substitution of 
the obtained value of a;* in Eq. (3) yields the corresponding 
real contact area A*. Alternatively, the value of a;* obtained 
from ECR measurements and Eq. (20) can be substituted in 
Eq. (2) to estimate the corresponding contact load and, in 
turn, the adhesion force and adhesion energy. This approach 
was adopted in a theoretical treatment of adhesion of homo- 
geneous conductive surfaces.*4 Because determination of the 
real contact area as a function of the contact load using Eqs. 
(2) and (3) involves surface topography parameters D and 
G*, the analysis is applicable for static contact conditions. 
Estimation of the real contact area using the previous ap- 
proach in the case of dynamic contacts can be accomplished 
by using the fractal parameters of the evolved (current) sur- 
face topography. This applies also to the alternative approach 
used to obtain the contact load as a function of the ECR from 
Eqs. (2), (20), and (21). 

Dividing Eq. (3) by Eq. (20) yields 


1/2 

Q AS 

ar=(2 ($5 exp(1.025ASg}”)(R*)7!. (22) 
L 

In view of Eq. (22), the real contact area can be directly 

obtained from simple ECR measurements that are indepen- 

dent of the surface topography and mechanical properties of 
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the contacting surfaces and the load applied, provided the 
drop in voltage is very small. The relationship between the 
real contact area and the ECR in the low-voltage regime [Eq. 
(22)] reveals dependence only on the properties of the thin 
insulating film. This is expected due to the dependence of the 
ECR on the microcontact area and film properties [Eq. (11)]. 
Although Eq. (22) was derived for either elastic or fully plas- 
tic microcontacts, it is applicable to the entire range of de- 
formation because of the intrinsic relationship between the 
contact area and the ECR existing at the single asperity level 
[Eq. (11)]. 

For the simple case of fully plastic microcontacts, an 
explicit relationship between the ECR and the contact load 
that is independent of the surface topography can be obtained 
for the low-voltage regime by using Eqs. (1) and (22), 


1/2 
P*= aia ($5 exp(1.025AS@;”)(R*)~!. (23) 
E \ pL t 

Equation (23) can be used to determine the adhesion energy 
and adhesion force in dynamic contacts, where the surface 
topography may change during operation due to deformation 
of the asperity microcontacts, similar to the method proposed 
for homogeneous conductive surfaces.’ 


IV. NUMERICAL ANALYSIS 


To obtain a general theory, it is necessary to account for 
elastic-plastic deformation of the asperity microcontacts”! 
and the voltage dependence of the tunnel resistance [Eqs. 
(12), (13), and (15)].” However, the complex constitutive 
relationships of the contact parameters in the elastic-plastic 
deformation regime and the nonohmic behavior of the tunnel 
resistance [Eqs. (12)—(15)], i.e., the dependence of the tunnel 
resistance of a microcontact on the voltage drop, inhibit the 
derivation of closed form solutions. Therefore, a numerical 
scheme”! was used to obtain results for the various contact 
parameters. The previous numerical scheme was updated in 
the present study to include an iterative procedure for the 
tunnel resistance. The procedure involves assuming an initial 
value of the voltage drop V across the contact interface [Fig. 
1(c)], and because all the microcontacts are connected in 
parallel and R,; is negligible compared to R,;, this is also the 
value of voltage drop V; across each microcontact due to the 
tunnel effect. The current Z; that flows through each micro- 
contact is then calculated according to the appropriate range 
of voltage [i.e., Eq. (8) or (12)]. Then, the total contact load 
and real contact area are obtained using Eqs. (5) and (6), 
respectively, and the total current J is calculated by 


N(ag) 


a I: (24) 
i=l 


Subsequently, the calculated total current is compared with 
the current applied and the assumed value of the drop in 
voltage is adjusted accordingly. This procedure is repeated 
until the variation in total current converges to a specified 
small tolerance value (e.g., 1%). Then, the total ECR is de- 
termined from the definition of the electrical resistance, 
R=V/I. 
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V. RESULTS AND DISCUSSION 


The numerical results presented in this section are for 
three-dimensional surfaces generated using fractal geometry 
for A,=1 wm’, 2<D<2.5,° 10° 7m<G<10 m (i.e., 
10° '<G*<10""), aj=1nm’, and 10°4<A*=<10-%, 
separated by a thin insulating film. The upper and lower 
limits of A* are selected so as to yield a small average size 
microcontacts relative to their average spacing in order for 
interaction effects to be secondary and extremely small real 
contact areas with unrealistically small number of asperity 
microcontacts not to occur. Thus, the curves shown in all the 
following figures are bounded by contact loads correspond- 
ing to A*= 107° and 1074. Unless otherwise stated, the re- 
sults presented below are for conductive surfaces with D 
=2.3, G*=10 7, E/Y=106, and v,=0.3, insulating film 
with t=10 A, K=6, and @j=2 eV, and applied current I 
=5 WA. 

Figure 2 shows the effect of fractal parameters and me- 
chanical properties on the variation of contact resistance R* 
with contact load P*. For given fractal parameters (i.e., 
fixed topography) and mechanical properties, the ECR de- 
creases with an increase in contact load due to the corre- 
sponding increase of the real contact area and the inverse 
dependence of the ECR on the real contact area [Eq. (22)]. 
Although for homogeneous conductive surfaces the ECR de- 
creases by about two orders of magnitude with an increase of 
the contact load,!” for conductive surfaces separated by an 
insulating thin film the ECR decreases by less than an order 
of magnitude over the same load range. This is due to inher- 
ent nonohmic behavior in the presence of an insulating thin 
film, where the tunnel resistance of each microcontact in- 
creases with a decrease in the drop in voltage across the 
microcontact.”’ Increasing the contact load results in lower 
ECR, which for fixed current flow J, yields a lower voltage 
drop. This, in turn, increases the microcontact tunnel resis- 
tance, thus compensating for the decrease in total ECR. This 
phenomenon is observed in the following figures showing 
the ECR as a function of the contact load and, hereafter, will 
be referred to as the voltage compensation mechanism. The 
fact that for a given contact load the ECR increases with the 
fractal roughness G* [Fig. 2(a)] is attributed to the smaller 
real contact area obtained with higher G* values (rougher 
surfaces). As discussed for isotropic conductive surfaces,” 
Fig. 2(a) can be used to determine G* from simple measure- 
ments of R*, provided all other parameters are known. Also, 
the ratio pet t can be determined from Fig. 2(a), Eq. (21), 
the given load, and the measured ECR. The value of P* 
determined from the measured ECR can be used to evaluate 
the adhesion force and interfacial adhesion energy following 
the method for homogeneous conductive surfaces described 
elsewhere.** The last two approaches can also be carried out 
by using any of the Figs. 2, 4, and 6 presented herein. 

Figure 2(b) shows the variation in ECR with the contact 
load and fractal dimension. For fixed load, increasing the 
fractal dimension causes the ECR to decrease. This is be- 
cause larger D values correspond to smoother (denser) sur- 
face topographies that produce larger real contact areas,!” 
which decreases the ECR. The almost identical ECR results 
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FIG. 2. Dimensionless electrical contact resistance R* vs dimensionless 
contact load P* for contacting rough surfaces: (a) G*=1077, 107°, and 
107!!, D=2.3, and E/Y=106, (b) D=2.1, 2.2, 2.3, and 2.4, G*=1077, 
and E/Y=106, and (c) E/Y=106, 288, and 391, G*=1077, and D=2.3 
(K=6, t=10A, gy=2 eV, and [=5 pA). 


obtained for D=2.1 and 2.2 are due to the very large critical 
truncated contact area a! [Eq. (4)] obtained with such small 
values of D, and result in fully plastic deformation of the 
majority of microcontacts. In this case the relationship be- 
tween the ECR and the contact load is independent of the 
surface topography [Eq. (23)]. Therefore, the performance of 
electrical contacts operating in the fully plastic deformation 
regime is insensitive to the surface topography, which may 
exhibit variations depending on the fabrication process 
and/or the evolution of wear. As for Fig. 2(a), the results 
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FIG. 3. Dimensionless contact area A* vs dimensionless electrical contact 
resistance R* for contacting rough surfaces with different mechanical prop- 
erties and fractal parameters (K=6, t=10 A, gj=2 eV, and J=5 uA). 


shown in Fig. 2(b) can be used to determine fractal dimen- 
sion D from the measured ECR and the contact load applied. 

Figure 2(c) shows the effect of the mechanical properties 
of the contacting surfaces on the variation of the ECR with 
contact load. For a given dimensionless contact load, in- 
creasing E/Y decreases the ECR. Since the reduced elastic 
modulus affects both E/Y and P*, it is difficult to draw a 
general conclusion about the effect of the mechanical prop- 
erties on the ECR. However, if the reduced elastic modulus 
is assumed to remain constant, the results shown in Fig. 2(c) 
can be used to obtain insight into the effect of the yield 
strength on the ECR. Because the mean contact pressure de- 
creases with the yield strength”! a larger contact area is re- 
quired to support a given contact load, resulting in lower 
ECR. Figure 2(c) can also be used to determine the yield 
strength of the softer material or the reduced elastic modulus 
(through an iterative procedure!” necessitated by the depen- 
dence of E/Y and P* on E) from ECR measurements. 

As discussed in Ref. 19 for homogeneous conductive 
surfaces, it is advantageous to obtain the real contact area in 
terms of the ECR because of the independence of this rela- 
tionship on the surface topography and mechanical proper- 
ties. This is also the case for conductive surfaces separated 
by a thin insulating film, as can be deduced from Eq. (22) for 
the low-voltage regime (V;=0). Figure 3 shows the real 
contact area as a function of the ECR for the entire range of 
fractal parameters and mechanical properties used to obtain 
the results shown in Fig. 2. The results were obtained by 
crossplotting data of R* (shown in Fig. 2) and A* (not 
shown here for brevity) obtained for the same contact load. It 
is interesting to note that all the data closely follow the same 
curve even for the general case of V;>0. This implies that 
the surface topography and mechanical properties do not in- 
fluence the dependence of A* on R*. This is of great impor- 
tance in dynamic systems where surface deformation may 
change the topography. 

According to the approximate analysis (Sec. IID), the re- 
lationship between A* and R* [Eq. (22)] is a function of t, 
K, and @ and is independent of the surface topography, 


L. Kogut and K. Komvopoulos 


mechanical properties, and contact load, in accord with the 
results shown in Fig. 3. However, significant differences ex- 
ist between the results obtained from Eq. (22) and those 
shown in Fig. 3 due to the assumption of V;=0 invoked in 
the approximate analysis. According to this analysis, the re- 
lationship between A* and R* is independent of the current 
used to obtain the ECR measurements, while obtaining an 
accurate solution depends on the current applied through the 
voltage drop that affects AS and z [Eqs. (13) and (15), 
respectively]. By curve fitting the numerical results shown in 
Fig. 3, the following empirical relationship was obtained: 


A*=3.21X 1077 exp(—2.91X 107 7R*). (25) 


The proportionality factor and exponent in Eq. (25) are both 
functions of t, K, @ 9, and the drop in voltage across the 
interface V, i.e., current J used in the ECR measurements. 

Figure 4 shows the effects of the dielectric constant, film 
thickness, and energy barrier height on the variation of the 
ECR with the contact load. For fixed contact load, the ECR 
decreases with K, t, and o, because R,; decreases with these 
parameters (Sec. IIC). For a given current, a smaller ECR 
produces a smaller drop in voltage across the interface, 
where the tunnel resistance is ohmic and, therefore, indepen- 
dent of the voltage drop. This may explain the pronounced 
variation of the ECR with the contact load in these cases, 
resembling the behavior of homogeneous conductive sur- 
faces with ohmic behavior.'? The results shown in Fig. 4 can 
be used to determine the dielectric constant, thickness, and 
energy barrier height of the insulating film from ECR mea- 
surements. However, because the dimensionless ECR is also 
a function of t and ọọ [Eq. (21)], an iterative procedure is 
required to determine these parameters. Hence, assuming an 
initial value of the unknown parameter (t or go), the dimen- 
sionless ECR can be calculated from the measured ECR us- 
ing Eq. (21). Then Fig. 4(b) or 4(c) can be used to determine 
whether the dimensionless ECR (based on the ECR measure- 
ments) coincides with the curve corresponding to the as- 
sumed value of the unknown parameter. This procedure can 
be repeated until eventually convergence to a specified toler- 
ance value is reached. 

The dependence of the real contact area on the ECR is 
shown in Fig. 5 for the ranges of dielectric constant, film 
thickness, energy barrier height, and contact load used to 
obtain the results presented in Fig. 4. The results were ob- 
tained using the same method as that used to obtain the re- 
sults shown in Fig. 3. As shown in Fig. 5, the relationship 
between A* and R* depends on the insulating film proper- 
ties. Equation (25) can be written in general form as 


A*=aexp(— BR*), (26) 


where œ and £ assume positive values that depend on K, t, 
o, and current J used in the ECR measurements. The values 
of œ and $ can be determined by curve fitting the numerical 
results obtained for the particular film properties shown in 
Fig. 5. Variation of the ECR with the real contact area is 
more pronounced for smaller values of K, t, and œo, due to 
the absence of the voltage compensation mechanism in these 
cases (see discussion of Figs. 2 and 4). Therefore, greater 
accuracy in the determination of the real contact area based 
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FIG. 4. Dimensionless electrical contact resistance R* vs dimensionless 
contact load P* for contacting rough surfaces: (a) K=4, 6, and 8, t 
=10A, and Yo=2 eV, (b) t=8, 10, and 12 A, K=6, and Qo=2 eV, and (c) 
o= 1, 2, and 3 eV, K=6, and t 10A (D=2.3, G*=107", E/Y = 106, and 
I=5 pA). 


on ECR measurements is expected for surfaces covered by 
relatively thin insulating films that possess both low dielec- 
tric constant and energy barrier. As for Fig. 3, the results 
shown in Fig. 5 are independent of the mechanical proper- 
ties, fractal parameters, and contact load. 

Figure 6 shows the effect of the current used in the ECR 
measurements on the variation of the ECR with the contact 
load. For fixed contact load, the ECR decreases with an in- 
crease in current because the tunnel resistance of a micro- 
contact decreases with an increase in the drop in voltage, as 
discussed in Sec. IIC. This implies that different applied 
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FIG. 5. Dimensionless contact area A* vs dimensionless electrical contact 
resistance R* for contacting rough surfaces with various values of K, t, and 
Po (D=2.3, G*=1077, E/Y=106, and I=5 uA). 


currents may bias the ECR measurements due to the voltage 
compensation mechanism. Another consequence of this ef- 
fect could be inaccurate estimation of the contact load. The 
results corresponding to 7=0O (also shown in Fig. 7) were 
obtained from Eq. (11). For current equal to 0.1 uA, the 
results obtained by using the voltage-dependent tunnel resis- 
tance relationships [Eqs. (12)—(15)] are in good agreement 
with the approximate results corresponding to /Æ=0. Rela- 
tively small differences occur at very light loads where the 
ECR assumes high values. In these light load cases, the high 
voltage drop invalidates the assumption of V;=0, even 
though the current is very low (J=0.1 uA). The results for 
I=0 and 0.1 uA were also found to be in good agreement 
with the approximate solution for fully plastic contacts [Eq. 
(23)], indicating the dominance of fully plastic microcontacts 
in lightly loaded contact interfaces. With such small currents, 
the ECR decreases by about two orders of magnitude due to 
the ohmic behavior of the tunnel resistance [Eq. (11)], while 
for higher currents the decrease in ECR is significantly less 
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FIG. 6. Dimensionless electrical contact resistance R* vs dimensionless 
contact load P* for contacting rough surfaces and various values of applied 
current J (D=2.3, G* = 10°’, E/Y=106, K=6, t=10 A, and @y=2 eV). 
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resistance R* for contacting rough surfaces with different mechanical prop- 
erties and fractal parameters, and various values of applied current J (K 
=6, t=10 A, and g)=2 eV). 


(for the same load range) due to the inherent dependence of 
the ECR on the voltage and the associated compensation 
mechanism. The results shown in Figs. 4 and 6 suggest that 
greater accuracy in the determination of the contact load 
from ECR measurements can be obtained for thin insulating 
films that possess both low dielectric constant and energy 
barrier and by applying a low current. 

The dependence of the real contact area on the ECR is 
shown in Fig. 7 for the entire range of electric current used to 
obtain the results presented in Fig. 6. The results were ob- 
tained with the same method used to obtain the results shown 
in Fig. 3. It should be noted that the relationship between A * 
and R* depends on the current used in the ECR measure- 
ments and follows Eq. (26). For current equal to or less than 
0.1 uA, the results for the voltage-dependent tunnel resis- 
tance [Eqs. (12)—(15)] are in good agreement with the ap- 
proximate results for the voltage-independent tunnel resis- 
tance, i.e., J=0 [Eq. (11)], except for high ECR values, in 
agreement with Fig. 6. The results corresponding to J=0 
(and 7=0.1 uA for low ECR values) are also in good agree- 
ment with the approximate solution for fully plastic contacts 
[Eq. (22)] (not shown here for clarity), confirming the domi- 
nance of fully plastic microcontacts in lightly loaded con- 
tacts. Figure 7 shows that it is advantageous to determine the 
real contact area from ECR measurements using low currents 
due to the more pronounced variation in ECR with real con- 
tact area at lower currents, a consequence of the ohmic be- 
havior of the ECR at low voltages. 


VI. CONCLUSIONS 


An ECR theory was introduced for conductive and rough 
(fractal) surfaces separated by a thin insulating film. The 
results illustrate the importance of an intervening film that 
consists of an insulating material, contact load, mechanical 
properties, surface topography, and current intensity in elec- 
trical contacts. The analysis accounts for current flow be- 
tween the surfaces by electron tunneling through a rectangu- 
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lar energy barrier that includes image forces. Based on the 
results and discussion presented, the following main conclu- 
sions can be drawn. 


(1) The constriction resistance plays a secondary role com- 
pared to the tunnel resistance and, therefore, the ECR is 
dominated by the tunnel effect. 

(2) The ECR decreases with an increase in contact load, 
fractal dimension D, and current flow, and a decrease in 
fractal roughness G, film thickness, dielectric constant, 
and energy barrier height. 

(3) Variation of the ECR with the contact load (or real con- 
tact area) is less pronounced than that observed in the 
absence of the insulating layer due to the intrinsic volt- 
age dependence of the tunnel resistance (nonohmic be- 
havior). For the same reason, this variation becomes less 
pronounced with an increase of applied current. 

(4) For contacting surfaces with small D values (i.e., D 
<2.2) the majority of the microcontacts are in the fully 
plastic deformation regime and the ECR is independent 
of the surface topography parameters. 

(5) For low currents (i.e., 7<0.1 uA), the numerical solu- 
tions for the ECR versus contact load and the real con- 
tact area versus ECR are in good agreement with closed- 
form analytical solutions obtained under the assumption 
of ohmic contacts and fully plastic deformation of the 
microcontacts. 

(6) The relationship between the real contact area and the 
ECR depends on the current, insulating film thickness, 
dielectric constant, and energy barrier height and is in- 
dependent of the contact load and the surface topography 
and mechanical properties of the contacting surfaces. 
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